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ABSTRACT 

It has long been expected that in some scenarios when a white dwarf (WD) grows to the Chandrasekhar 
limit, it can undergo an accretion induced collapse (AIC) to form a rapidly rotating neutron star. Nevertheless, 
the detection of such events has so far evaded discovery, likely because the optical, supernova-like emission is 
expected to be dim and short-lived. Here we propose a novel signature of AIC: a transient radio source lasting 
for a few months. Rapid rotation along with flux freezing and dynamo action can grow the WD's magnetic 
field to magnetar strengths during collapse. The spindown of this newly born magnetar generates a pulsar wind 
nebula (PWN) within the ~ 10~ 3 - 10 _1 M Q of ejecta surrounding it. Our calculations show that synchrotron 
emission from the PWN may be detectable in the radio, even if the magnetar has a rather modest magnetic 
field of - 2 x 10 14 G and an initial spin period of ~ 10ms. An all-sky survey with a detection limit of 1 mJy at 
1 .4 GHz would see ~ 4(//10~ 2 ) above threshold at any given time, where / is the ratio of the AIC rate to Type 
la supernova rate. A similar scenario may result from binary neutron stars if some mergers produce massive 
neutron stars rather than black holes. We conclude with a discussion of the detectability of these types of radio 
sources in an era of facilities with high mapping speeds. 

Subject headings: stars: magnetic fields — stars: neutron — stars: winds, outflows — white dwarfs 



1. INTRODUCTION 

As an accreting white dwarf (WD) grows toward the 
Chandrasekhar limit, a well-known potential outcome is 
ignition of its n uclear fuel, leading to a Typ e la su- 
pernova (SN la, iHillebrand t & Niemeyer 2000). How- 
ever in some cases (e.g. mass transfer onto O/Ne/Mg 
WDs and C/O WD m ergers; ICanal & Schatzmanl 1 19761: 
iNomoto & Kondolll991l) electron capture can rob the core 
of its degeneracy pressure support leading to formation of 
a neutron star (NS). This "Accretion Induced Collapse" 
(AIC) has been invoked to explain m illisecond pulsars (e.g. 
Bhattacharya & van den Heuvelll99l|). subsets of gamma-ray 
burst s (e.g. Par et al.lll992t iMetzger et al.ll2008bl) . magnetars 
(e.g. IUsovl ll99ir and may be a source of r-pro cess nucle- 
osynthesis (IHartmann eUull985tlFrver et al.lll999l) . 

Despite its potential importance, there has been no reported 
detection of an AIC event. To start with, th e expected AIC rate 
is no more than » 1 % of that of SNe la ( Yunge lson & Liviol 
119981) . Next, relative to Type I and Type II SNe, the ejecta 
mass is expected to be small (< 10 _1 M Q ), produce little 56 Ni 
(< 10~ 2 M Q ), and move at high velocity (w 0.1c). The result- 
ing optical transient is thus considerably fainte r than a typical 
SN (5 magnitudes or more ) and lasts ~ 1 day (IMetzger et al.1 
2009bl iDarbha et al.ll2010h . 

Rapid rotation should accompany AIC due to the accre- 
tion of mass and angular momentum. Furthermore, a strong 
magnetic field may be amplifie d through flux freezing durin g 
collapse and via dynamo action dDuncan & Thom psonl 19921) . 
Therefore a plausible outcome of AIC is creation of a quickly 
spinn ing magnetar (lUsovll 19921: Iking et al.ll200ll iLevan et al. 1 
2006). The WD collapse unbinds material (iDessart et alJ 



2006), and the remnant disk loses mass via outflows driven 
by neutrino heating, turbulent viscosity, and recombina- 
tion of free nuclei into helium dL"ee & Ramirez-Ruizl 120071 : 
IMetzger et all l2008al I2009at iLee et all 120091) . This leads to 
M e j » 10~ 3 - 10 _1 M Q of ejecta with velocity v e j ~ 0.1c. Such 
a configuration will also follow a double NS merger if the to- 



tal binary mass is below the maximum NS mass. 

Here we investigate the detectability of the pulsar wind neb- 
ula (PWN) powered by a newly formed magnetar following 
AIC. In $2] we describe a model for how the PWN expands 
into the surrounding ejecta. In $3] we estimate the radio spec- 
trum and lightcurve, showing that the PWN is a transient radio 
source for a few months. We estimate the detection rate in $4] 
and discuss the detectability of such events with soon-to-be- 
commissioned high speed radio mapping machines in §|5] 

2. DYNAMICS OF THE PULSAR WIND NEBULA 

Following the implosion of the WD and subsequent disk 
outflows, the ejecta expands with velocity v e j « 0.1c and 



kinetic energy E e j 



M ej v|/2 



10 erg. This plows into 



the ISM with particle density no, which can vary greatly 
from ~ 10~ 6 cm" 3 for events outside of their host galaxy to 
~ 1 cm" 3 in denser environments. Once the ejecta has swept 
up a mass comparable to its own, it de celerate as it enters 
the S edov-Taylor phase on a timescale dMcKee & Truelovel 
[19951) 



:0.5£ e -: 1/2 M e 5 , /6 K"or 1/3 



1/CZ7- 1 / 2 j.f5/6 -1/2 

16E 50 < MJ 2 n yr, 



(1) 

where £50 = £ej/10 50 erg and M_ 2 = M ej /1O" 2 M . Since r ST 
is much longer than the times of interest (especially when «o 
is small), the ejecta is always in an ejecta-dominated (or free- 
expansion) phase. In this case the ejecta maintains constant 
velocity with radius « v e jf . The forward shock is merely a 
distance « (4-/3)^ 3 R e j ahead of R e j, and the reverse shock has 
barely developed behind it. The pressure behind the forward 
shock and down to the reverse shock is roughly constant and 
given by the strong shock limit w (A/3)m p nov\y 

The magnetar injects energy in the form of magnetic fields 
and relativistic particles at a rate 



L(t) = L i) /(l+t/T)P 



(2) 
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Here we assume dipole spindowrQ and thus p = 2. For a mag- 
netic moment \i, initial spin frequency fio = 2tv/Pq, and mo- 
ment of inertia / = 0.35M*/? 2 , dLattimer & Pr akash 200j), 



L = fi 2 n 4 /6c 3 : 



1.2 x lO^ls^ergs" 



(3) 



and 



r = 6Ic s /(/Sll = 4.8 x IOV33P3 s 



2n2, 



(4) 

where /133 = /i/10 33 G cm 3 , F3 = Po/3ms, and we use M* = 
1.4M Q and/?* = 12km. 

The energy input from the magnetar powers a PWN with 
radius R p . Deeper inside, a wind termination shock forms at 
radius R,, where the ram pressure of the pul sar wind equals 
the PWN pressure P (Gaensler & Slane 2006b, 



R r » (L/4ircP) 



1/2 



(5) 



where R t <C R p - The general picture developed from galactic 
PWNe is that particles are accelerated at or near R, . This seeds 
the PWN with relativistic electrons out to R p . 

A central issue is the content of the PWN. As a pulsar wind 
flows from the light cylinder, it is inferred to have a large mag- 
netization (cr ~ 10 4 , where a is t he ratio of Poynting flux to 
particle energy flux; Arons 2002). However multiple lines of 
evidence, including the expansion velocities and high-energy 
modeling of PWNe, require that the magnetization decreases 
significantly (cr 1) by the time the wind reaches R,. In the 
region where the radiation arises, R, < r < R p , we define r\ e 
and r\B as the fraction of the magnetar luminosity that goes 
into electrons and magnetic fields, respectively, with typical 
values of r\ e ~ 0.999 and t)b ~ 10~ 3 . 

Following lRevnolds & Chevalier! (119841) we assume the en- 
ergy density of electrons and magnetic fields evolve indepen- 
dently and obey a relativistic equation of state. Thus 



j t (4nP e Rl) = (r ]e L-A)R p , 
where A is the radiative loss rate, and 



VuLRp, 



(6) 



(7) 



where Pg = B 2 /%tt. Note that in the adiabatic limit r/^L w 0, 
this predicts Pb oc R~ p and B oc R~ 2 , as expected from flux 
freezing. Momentum conservation is given by 

P-pJ Vp -^L )} , (8) 



M S ^=A,R] 



where v p = dR p /dt and P = P e +P B is the total pressure. 
The PWN sweeps up ejecta and creates a shell with mass 



M s ■ 



M ei (R p /v ei t) 3 , R p < v ej f 



Rp > V e jf ' 



(9) 



Once R p « v e jf, the PWN reaches material that has been 
shock-heated by its interaction with the ISM. Since we gen- 
erally find P ^ (4/3)m / ,noi , ej' we d° not ex P ect the shock- 
heated pressure to have a large dynamical effect on the PWN. 
This is in contrast to PWNe growing within a SN rem- 
nant, where the large pressure behind the reverse shock (on 

1 When the magnetar is spinning Pq < 3 ms and readily radiating neu- 
trinos, this can give r ise to a neutrino-driven, magneto-centrifugal wind 
I Thompson et al. 2004) which could greatly enhance the spindown. 



a timescale t ^> tst) can ca use significant compression and 
magnetic field amplification (Reynolds & Chevalier 1984). 

Equations (O, (IT), and (O with dR p /dt = v p provide four 
first-order differential equations for the four dependent vari- 
ables P e , P B , R p , and v p , respectively. For simplicity, we cal- 
culate the PWN evolution using A 0, so that equations © 
and (0 are combined using P = P e +P B and rie + VB = 1 ■ When 
L is constant, the analytic result is (Chevalier 1977b 



R p « 1.3 x 10 16 L;f4 /I0 M:^ 2 f6 6/5 cm, 



and 



1/5 r l/5 P -9/20 M 3/4 t -13/10 



50 



2 '6 



G, 



(10) 



(11) 



where L 47 = L/10 47 erg s _1 , f 6 = f/10 6 s, and r] B .-i = Vb/10~\ 
This gives some idea of the rough values expected, although 
for our detailed calculations, R p and B deviate slightly from 
these scalings when t > r. 

Rayleigh-Taylor instabilities can act at R p due to the low 
density PWN pushing up against high density ejecta. The 
growth time for a large density contrast is trj ~ (g e sk)~ 1 ' 2 , 
where g e s is the effective gravitational acceleration at the 
boundary and k the wavenumber. For R p given by equa- 
tion ([Toll, geff ~ d 2 R p /dt 2 w (6/25)R p /t 2 and the growth rate 
for k « n/Rp is trt ~ (25/6n) I//2 f. For sufficiently small 
wavelengths (large «), trt < t and instability results. This 
is not surprising since similar systems, like the Crab Nebula, 
have morphologies strongly impacted by instabilities. But for 
the current analysis we do not include this complication. 

3. SYNCHROTRON RADIO EMISSION 

A power-law spectrum of relativistic electrons n(E) = KE~ S 
are accelerated near the termination shock at R, and fill the 
PWN out to R p , where n{E) is in units of electrons erg" 1 cm -3 . 
The power-law parameters within the PWN c an change via 
cooling and injection of new electrons (e.g. Gelfa nd et al.l 
120091: Buccianti ni et al.|[201 ll) . For the present work we es- 
timate the synchrotron spectrum at any time by fixing P e and 
Pb from our dynamical calculations. 

Our d iscussion of synchr otron emission largely follows the 
work of Pacholczyk ( 1970). Synchrotron emission is self ab- 
sorbed below a frequency 

zaja = 2 Cl (/? ;)C6 ) 2 / ( ' 5+4) J rr 2 / ( ' 5+4) (Bsin0) (l+2) / (j+4) , (12) 

where c\ = 6.27 x 10 18 i n cgs units, cj, depe nds on s and can be 
found in Appendix 2 of Pacho lczvkl d 19701) . and 8 is the pitch 
angle. Throughout we assume an average of sin 9 = (2/3) 1 ' 2 . 
In the optically thick limit the flux is 

5/2 



^£5 

D 2 c 6 



(Bsinfl)- 1 / 2 ( — 
2c 1 



(13) 



where D is the distance and C5 is another constant. In the 
optically thin limit, 

-0-i)/2 



F, = -^c 5 K{Bsmef S+1 V 2 

3D £ \ 2ci 



(14) 



which assumes an emission filling factor of order unity. We 
use a simple interpolation between these two limits for the 
total emission spectrum F v . The location of the spectrum's 
peak at ^sa oc. Z}( j+2 V( s + 4 ) shifts to lower frequencies as B de- 
creases during expansion (eq. ifTTI ). The peak flux scales 
oc b( 2s+3 )A i+4 \ and thus also decreases with time. 
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The synchrotron cooling time scales as E/\E\ oc u~ l l 2 . 
Therefore there is a v c above which synchrotron c ooling beats 
adiabatic expansion (Reynolds & Chevalier 1984) 



c 2 fi 3 sin 3 9 \R 



(15) 



where c-i = 2.37 x 10~ 3 in cgs units. Electrons that evolve adi- 
abatically maintain their spectrum, while electrons that cool 
from synchrotron emission steepen, therefore 



«(£): 



KE~ S , E < E c 
KE c E~ (s+l \ E>E C > 



(16) 



where E c is the energy of electrons emitting a frequency v c . 
The prefactor K is set by the total energy density of electrons 
U e = 3P e , found by integrating 



U e = I n(E)EdE, 



(17) 



where E m i n and E m!a are the minimum and maximum ener- 
gies of electron spectrum, respectively. For s ~ 1.5 and the 
limit E m - dx 3> E c ^> E^, U e ~ 4KE 2 ~ S . The energy density is 
roughly set by E c since the spectrum steepens for E > E c . 

We next consider p otential free-free absorptio n. The ab- 
sorption coefficient is (Rvbicki & Lightman 1979) 



1.9 x 10" 2 r" 3/2 Z 2 n e n,z/" 2 g ff cm" 



(18) 



for hv <C k B T, where T is the temperature, Z is the average 
charge per ion, n e and n -, are the electron and ion densities, re- 
spectively, gff ~ 1 is the Gaunt factor, and all quantities are in 
cgs units. Just behind the forward shocks (see the temper- 
ature is high (3> 10 9 K), and free-free absorption is negligible. 
But at the front edge of the ejecta, with density p e j and tem- 
perature T, pressure continuity across the contact discontinu- 
ity requires p e jk B T jm p ~ notn p v^. From this we estimate that 
T ~ 10 5 - 10 8 K at t ~ 10 6 - 10 7 s. For Z/A « 1 /2 (where A is 
the average atomic number) we estimate Z 2 « e «, ~ (p e j/m p ) 2 . 
The observed flux is thus F u o \, s w F u e~ T ", where Tff w AR p 
and AR p w 0AR p is the shell thickness. 

In the top panel of Figure [TJ we plot the evolution of the 
key frequencies v^a (eq- lfl2"l ). v c (eq. |[l"5l ). and ug, where 
the latter is defined by tq ■ « 1 . From this one can follow the 
spectrum peak (at i>sa) and determine when it is detectable. 
For example, 10GHz emission peaks at w 6 x 10 6 s. In the 
bottom panel we explore the lightcurves at 1 .4 GHz as we vary 
T and t]b- Typical timescales are ~ months with a peak of 
~ 3- lOmJy. Free-free absorption is mostly negligible unless 
T < 10 4 K. The peak flux and time of peak are sensitive to t]b- 

In Figure |2] we plot the peak flux at 1.4 GHz for a range 
of dipole field strengths and initial spin periods. The general 
trend is that faster spins and larger fields result in a more lumi- 
nous radio source. This reverses in the top left corner where 
r is sufficiently short that L(t) has decreased significantly by 
the time of peak. Over most of this parameter space, the mag- 
netar winds are not e xpected to be strong eno ugh to generate 
a collimated outflow (Bucciantini et al. 20 12|). 

4. DETECTION RATES 

AIC can result from both channels popularly discussed for 
SNe la (single-degenerate and double-degenerate; see ®. 
Thus it makes sense to normalize the AIC rate to that of the 
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FIG. 1. — The top panel shows the evolution of the critical frequencies 
with time for M e j = 1CT 2 Mq, v c j = 0.1c, fi = 10 33 Gcm 3 , and P = 3 ms. 
The electron spectrum has a power-law index s = 1.5. The shaded region 
shows where free-free absorption suppresses the synchrotron spectrum for 
T = 10 5 K. This moves up for lower T (see eq. 1181 ). The spectrum peaks 
at frequency i/sa- The bottom panel shows the time-dependent flux at v — 
1.4GHz for a distance D = 100 Mpc over a range of parameters as labeled. 
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FIG. 2. — The peak flux at 1.4 GHz at 100 Mpc as a function of the magnetar 
field B t and its initial period Pq. The contours of constant peak flux are 
spaced logarithmically in units of mjy with the thick, solid curves marking 
the labeled contours. In all cases we use r)g = 10~ 3 , T = 10 5 K, and * = 1.5. 

SN la rate. The Lick Observatory Superno va Search find s 
a rate of (3.01 ± 0.062) x 10" 5 Ia Mpc" 3 yr" 1 dLi et alj|201 ll 
which corresponds to (4.0-7.1) x 10~ 3 Ia yr" 1 for the Milky 
Way. Using population synthesis, Yungelson & Livio ( 1998) 
find AIC rates of 8 x 10" 7 - 8 x 10" 5 AIC yr" 1 for the Milky 
Way, depending on assumptions about the common-envelope 
phase and mass transfer. Their upper bound is similar to the 
constrain t obtained from observed abundances of neutron rich 
isotopes dHartmann et al.|[l985l : lFrver et al.| [l999). 
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For a typical peak 1.4GHz flux of f^ioo ~ 5mJy at a dis- 
tance of lOOMpc, and an AIC rate that is a fraction / of the 
SN la rate, the detection rate of radio transients from AIC is 

^<>^™>~KT^)(i§r 2yr ~' <i9) 

where / ~ 10~ 4 - 10~ 2 . For a ~ 3 month duration emitting 
above ~ 1 mJy, we expect ~ 4(//10~ 2 ) AICs above threshold 
at a given time. In contrast, merely a few AICs are estimate d 
to be detected as kilonovae per year (lMetzgeretal.ll2009bl) . 
and this number could be muc h less if a l ack of differential 
rotation during WD a ccretion dPirol 12008)) in hibits disk for- 
mation upon collapse (Abdikamalov et al. 20101) 

Similar radio emission is possible if some NS mergers 
produce a massive NS. This would then provide an electro- 
magnetic counterpart following mon ths after the gravitational 
wave emission durin g coalescencf] (Metzger & Berger fl2bl2t 
iNissanke et al.l l2012). In fact, the distance at which the radio 
emission can be detected is similar to that probed by the next 
generation of "advanced" ground-based laser-interferometers. 
The NS merger ra te is comparable or greater than the AIC 
rate. For example, iKim et afl (120051) estimate a Galactic rate 
of ~(0.1-3)x 10~ 4 yr _1 . The detection frequency is thus sim- 
ilar to equation (fT9l with / = 10~ 2 , with the caveat that it may 
be much less if most NS mergers produce black holes. 

5. DETECTION & IDENTIFICATION OF AIC EVENTS 

This is a timely topic to discuss given the renaissance that 
is now occurring in decimetric radio astronomy. Refurbished 
(EVLA) and new facilities promise high mapping speeds by 
using small diameter antennas to realize a given total area (e.g. 
MeerKAT), through the use of focal plane arrays instead of a 
single feed to gain massive multiplex advantage (APERTIF), 
or both (ASKAP). 

We offer a plausible project aimed at detection of AICs that 
can be undertaken over the next few years. For a la detec- 
tiorQ of a 1 mJy point source, the EVLA mapping speed in 
the 1 .4 GHz band (for a 500 MHz band width) is ss 86?? square 
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